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2 ■ ABSTRACT 

> ■ 

We studied the metal-poor globular cluster populations of a large variety of galaxies 
J^j ' and compared their mean metallicity with the properties of the host galaxies. For this 

purpose, we constructed a comprehensive database of old, metal-poor globular cluster 
populations, hosted by 47 galaxies spanning about 10 magnitudes in absolute bright- 
ness. The mean metallicities of the systems are found to be very similar and to lie 
in the —1.65 <[Fe/H]< —1.20 range (74 % of the population). Using only globular 
cluster systems with more than 6 objects detected, we find that 85 % of the popula- 
tion are within —1.65 <[Fe/H]< —1.20. The relation between the mean metallicity of 
the metal-poor globular cluster systems and the absolute V magnitude of their host 
galaxies presents a very low slope which includes zero. An analysis of the correlation of 
the mean metallicity of the populations with other galaxy properties (such as velocity 
dispersion, metallicity, environment density) also leads to the conclusion that no strong 
correlation exists. The lack of correlation with galaxy properties suggests a formation 
of all metal-poor globular clusters in very similar gas fragments. A weak correlation (to 
be confirmed) might exist between mean metallicity of the metal-poor clusters and host 
galaxy metallicity. This would imply that at least some fragments in which metal-poor 
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globular clusters formed were already embedded in the larger dark matter halo of the 
final galaxy (as oppose to being independent satellites that were accreted later). Our 
result suggests a homogeneous formation of metal-poor globular clusters in all galax- 
ies, in typical fragments of masses around 10 9 -10 10 M with very similar metallicities, 
compatible with hierarchical formation scenarios for galaxies. 

We further compare the mean metallicities of the metal-poor globular cluster pop- 
ulations with the typical metallicities of high redshift objects. If we add the constraint 
that globular clusters need a high column density of gas to form, Damped Lyman a sys- 
tems are the most likely sites among the known high redshift objects for the formation 
of metal-poor globular cluster populations. 

Subject headings: Galaxies: formation; Galaxies: halos; Galaxies: star clusters; Cosmol- 
ogy: observations; Cosmology: early Universe 

1. Introduction 

Globular clusters contain the oldest known stellar populations of the Milky Way and proba- 
bly of the observed Universe. Consequently, they hold cosmologically significant information and 
are often used as fossil records of the formation of galaxies. At least two globular cluster sub- 
systems were found to coexist in our Galaxy (Kinman 1959; Zinn 1985; Armandroff & Zinn 1988; 
Armandroff 1989): i) a metal-poor, slowly rotating spherically distributed one in the halo and 
ii) a metal-rich, rapidly rotating population concentrated in the disk (or the bulge, e.g. Minniti 
1995). A similar situation was observed in M 31 (Ashman k, Bird 1993; Barmby et al. 2000). 
Subsequently, globular cluster su6-populations were also discovered in early-type galaxies (Zepf & 
Ashman 1993). This motivated comparisons between the globular cluster sub-populations and the 
host galaxies. Such studies were pioneered by Ashman & Bird (1993) and revived by Forbes et 
al. (1997a). Our work is motivated by the recent increase in data identifying metal-poor globular 
cluster populations, especially in early-type galaxies, as well as by recent studies associating the 
metal-poor sub-population in early-type galaxies with a extended, halo-like(?) component, similar 
to the ones observed in late-type galaxies (Geisler et al. 1996; Kissler-Patig et al. 1997a; Lee et al. 
1998). This encouraged us to look not only into common properties of these sub-systems and to 
investigate possible correlations between their properties and those of their host galaxies, but also 
to put them into a galaxy formation context and to identify their possible sites of formation. 
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2. The compilation of old globular cluster populations 

2.1. Focusing on metal-poor globular clusters 

We intend to select the oldest globular clusters around galaxies, with the minimum pollution 
from younger globular clusters. Even if absolute ages are not well defined, halo globular clusters in 
the Milky Way were shown to be very old systems and at least older than 10 Gyr (e.g. Chaboyer et 
al. 1998; Jimenez &: Padoan 1998). Globular clusters with [Fe/H] < —1.2 are essentially located in 
the halo but, more importantly, they are likely to be coeval within less than 1 Gyr (see Rosenberg 
et al. 1999). This also seems to be the case in other Local Group galaxies (Olszewski et al. 1996; 
Sarajedini et al. 1998). Furthermore, spectroscopic studies of giant ellipticals showed that the age 
of their metal-poor globular clusters appears indistinguishable within the errors (1 to a few Gyr) 
from that of the Milky Way halo clusters (Kissler-Patig et al. 1998a; Cohen et al. 1998). Assuming 
that these globular clusters are among the first stellar populations formed in the galaxies, we expect 
them to reflect the local pre- or proto-galactic conditions and especially the abundances more than 
10 Gyr ago. This paper will focus on the metal-poor globular clusters, assumed to be the oldest 
ones. 

It was proposed that, in some galaxies, the metal-poor and (at least part of) the metal-rich 
clusters were coeval within a few Gyr (e.g. Ortolani et al. 1995; Feltzing & Gilmore 2000; Puzia et 
al. 1999). However, given the fact that this has not been generally demonstrated, and given that 
a number of scenarios predict the metal-rich clusters to be younger than the metal-poor ones, we 
will not discuss these in this paper. 



2.2. On colors and metallicities and the detection of bimodality 

Our goal is to identify properties that are commonly measured in all known metal-poor pop- 
ulations. The mean metallicity of the metal-poor sub-population in a galaxy (noted as [Fe/H] mp 
throughout the paper) is the only such property currently available for a majority of galaxies. The 
following analysis will focus on this property, and by mean metallicity we actually refer to the peak 
of the metallicity distribution of the metal-poor globular clusters sub-population in a galaxy. 

The determination of the mean metallicity of metal-poor populations is complicated by two 
problems when derived from broad-band colors, as for the majority of our sample: i) a perfect 
separation of the metal-poor and metal-rich populations (and thus a determination of the mean 
color /metallicity for an unbiased sample of metal-poor clusters) is not possible from the broad- 
band color distributions, and ii) the transformation of the broad-band colors into metallicities 
suffers from systematic errors. The uncertainties associated with each point are discussed below. 

In order to determine the mean metallicity of the metal-poor clusters from a color distribution 
of a whole system, the distribution is probed by a mixture modeling test (KMM, see Ashman et 
al. 1994) that, among others, returns the most likely color peaks of the sub-populations. Typical 
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errors for the peak determination in (V — I) (the most widely used color), induced by the KMM 
method alone, range from 0.01 to 0.05 magnitudes (sensitively dependent on the number of data 
points and the intrinsic width of the distributions) which translates into errors in the mean [Fe/H] 
values of up to 0.25 dex. These errors are typically added to the statistical errors present in the 
photometry, and to potential systematic errors in the sampling of the metal-poor population. A 
large fraction of our metallicities are obtained from (V — I) colors derived from WFPC2/HST data. 
Median color values for the same system can vary from author to author by up to 0.06 in (V — I) 
(see for instance NGC 4472 in Neilsen & Tsvetanov 1999; Puzia et al. 1999). Moreover, in order 
to combine the studies in different bands, and to combine the results derived from colors and from 
spectroscopy, the broad-band colors need to be transformed into metallicities. The sensitivity of a 
color to metallicity transformation varies by almost a factor of two when going from (V — I) over 
(B — I) to (C — Ti), making a homogeneous compilation difficult. The different transformations for 
a given color into metallicity (often derived from the Milky Way clusters) can introduce errors of 
the order of [Fe/H]~ 0.3 dex depending on the exact method used to derive the relation (e.g. the 
comparison in Kissler-Patig et al. 1998b). 

In summary, taking the quoted mean metallicities directly from the literature could result in 
a artificial scatter of up to 0.4 dex in [Fe/H] mp in the extreme cases, given the different analyses of 
the various authors. Therefore, in addition to the [Fe/H] mp values, we will use the (V — I) values 
of the blue peak in the globular cluster color distributions, available for a subset of our sample. 
The latter avoids the possibility of introducing any error related to a different method of deriving 
metallicities, or errors associated with the conversion of colors to metallicities. 

2.3. The data compilation 

The mean metallicity of the metal-poor clusters, as well as their mean color (when available) 
are given in Table 1 for all galaxies which are known, to date, to host a distinct metal-poor cluster 
population. We included the data published by Kundu in his PhD thesis (Kundu 1999) for 9 new 
galaxies (and additional data for 8 galaxies already in our list). 

New values for 7 galaxies, derived from data presented in Paper I (Gebhardt & Kissler-Patig 
1999), are also added. The data and reduction are presented in the original paper. We selected 
all galaxies with clear bimodalities (see Paper I), and used the KMM code to derive the peak 
color of the blue globular clusters. These values, as well as the ones from Kundu (1999) were 
transformed into metallicities using the ([Fe/H], (V— I) relation given in Kissler-Patig et al. (1998b): 
[Fe/H]= 3.27 ■ (V — I) — 4.50. 

All other values in Table 1 are taken from the original references. We also added to the above 
sample a number of dwarf galaxies which present a unimodal metallicity distribution function 
with an average metallicity below the threshold defined in Sect. 2.1. We have assigned global 
uncertainties to each color : a(y_j^ ~ 0.25, &(b-i) ~ 0.20, g^c-t-i) ~ 0-15 (following Forbes et al. 
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1997a) unless different values are given in the original papers. 

The compilation includes galaxies of all types. However spiral galaxies are under-represented 
and bright elliptical galaxies dominate the sample. This observational bias is essentially due to the 
fact that : i) globular clusters are more easily identified on the smooth background of elliptical 
galaxies than in dusty spirals, it) bright ellipticals host a larger number of globular clusters than 
faint ellipticals or spirals. 

The host galaxy properties to which we compare the globular cluster properties are compiled in 
Table 2. These are taken from the HYPERCAT database (Prugniel & Simien 1996; Golev & Prugniel 
1998), except for the environment density taken from the Nearby Galaxy Catalog (Tully 1988). We 
chose the Mg2 index as a metallicity indicator rather than, for example, the color of the galaxy, 
because it is widely available and does not require dereddening or transformation to reflect the 
metal content of the galaxy. We keep in mind that Mg2 is essentially measured in the very inner 
regions of the galaxy, and does not directly reflect the mean metallicity of the halo. Nevertheless, it 
appears to be a good indicator for the final global metallicity and correlates well with the velocity 
dispersion of the galaxy (e.g. Dressier et al. 1987). 

The above mentioned velocity dispersion is used as a size indicator for the galaxies, rather than 
e.g. the estimated absolute magnitude, since the former is distance independent, while the latter is 
not and is therefore more difficult to bring onto a homogeneous scale. In order to get comparable 
values for all galaxies, we tried to select only seeing-limited ground-based determinations for the 
central velocity dispersion, e.g. HST/STIS values being systematically higher due to the higher 
spatial resolution. 

The data described above is used in the next section to investigate possible correlations between 
the mean metallicity of the globular clusters and the host galaxy properties. 

3. Globular cluster mean metallicities and galaxy properties 

3.1. A "universal" mean metallicity for metal-poor globular clusters 

Until the early 90s, the mean metallicities of the globular cluster systems was thought to 
correlate with the galaxy luminosity (vandenBergh 1975; Brodie & Huchra 1991). Ashman k, Bird 
(1993) first investigated the correlation between the mean metallicity of the metal-poor globular 
clusters only and the galaxy luminosity. They based their analysis on data of 4 local dwarf galaxies 
as well as the LMC, the Milky Way and M31. They found that a mean value of [Fe/H]~ —1.6 dex 
for all halo globular cluster systems appeared consistent with the data, and claimed that genuine 
halo globular cluster systems have comparable mean metallicities, irrespective of the parent galaxy 
luminosity. Further, adding the data of the 4 early-type galaxies known at that time to have bimodal 
globular cluster metallicity distributions, they speculated that the earlier relations (see above) were 
primarily a result of the high fraction of metal-rich globulars in bright elliptical galaxies. 
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Forbes et al. (1997a) later confirmed this result with a slightly larger sample of 11 galaxies, by 
looking at the correlation of mean metallicity with galaxy luminosity for the metal-poor and metal- 
rich clusters separately. For the metal-rich population, they found a positive correlation at the 3 a 
level. For the metal-poor globular clusters, Forbes et al. (1997a) did not detect any correlation but 
a random scattering. Their dataset has a mean of [Fe/H] mp = —1.16 with a dispersion of a = 0.28 
(as determined by us using a maximum- likelihood estimator on their 11 values). 

The sample of globular cluster systems presented in this paper is the largest database to-date, 
and > 4 times more numerous than Forbes et al. (1997a) initial dataset. The mean of our data lies at 
[Fe/H] mp = —1.45 with a dispersion of a = 0.15. Compared to the Forbes et al. dataset, our sample 
is slightly more metal-poor on average and exhibits a smaller scatter. The mean absolute magnitude 
of the galaxies in our sample is < My >= —20.1 with a dispersion of 2.4 mag. Fig. 1 shows the 
distribution of mean metallicities, plotted as a percentage of globular cluster systems within each 
bin (A[Fe/H] mp = 0.15). The first apparent result is that the mean-metallicities of the metal-poor 
globular clusters are not distributed homogeneously over the spanned range, but rather peak around 
a characteristic value : 74% of the sample is concentrated around —1.65 <[Fe/H] mp < —1.20. The 
distribution is asymmetric and more extended towards lower metallicity (dominated by the dwarf 
galaxies), resembling in shape to the halo field star metallicity distribution. When we exclude 
galaxies with a number of observed blue globular clusters Nqc < 6 (almost exclusively dwarf 
galaxies) from the sample, our statistics are slightly altered to [Fe/H] mp ~ —1.37 ± 0.2 and for the 
39 remaining galaxies < My >= —20.8 ± 1.6. The new histogram (Fig. 1) appears more peaked. 
Indeed, 85 % of the globular cluster systems are found within —1.65 <[Fe/H] mp < —1.20 and 67% 
within the two central bins i.e. —1.55 <[Fe/H] mp < —1.25. Thus, the data suggest an even higher 
concentration of the metal-poor globular clusters. 

The peaked (roughly Gaussian) distribution of the mean metallicities would be expected on 
the base of the Central Limit Theorem if all metal-poor globular clusters could be associated with 
a single sample (i.e. be considered to have a similar origin). The fact that the distribution looks 
indeed peaked supports the latter hypothesis. 

Figure 2 shows the metallicities of the globular cluster systems as a function of the abso- 
lute magnitude My. The global trend is a decrease of the metallicity with My. The statis- 
tical Spearman's rank test seems to confirm this impression and gives a probability of 0.0005 
that a correlation is not present (Spearman's p = -0.521). A linear fit gives : [Fe/H] mp = 
— 0.06(±0.01)My — 2.72(±0.22) for 46 values. However, removing the globular cluster systems 
with Ncc 1 < 6 from the sample (but keeping Kundu's), the same Spearman rank test gives a much 
lower probability of 0.0813 that a correlation is not present (Spearman's p = -0.283 for 39 values). 
A linear fit gives [Fe/H] mp = -0.02(±0.02)M F - 1.87(±0.31), with a low slope not significantly 
different from zero (at la). 

Thus, the mean metallicity of the old, metal-poor globular clusters seems to correlate with 
the absolute luminosity of their host galaxy. However, taking only the galaxies with N^c > 6, 
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this correlation is no longer statistically significant while we still have a large range in galaxy 
luminosity ( — 23 < My < —16). Thus, our findings confirm Ashman $ Bird's results for the metal- 
poor clusters. However, the globular cluster systems of dwarf galaxies seems to deserve a more 
complete discussion and they will be discussed in a future paper. 

3.2. Mean metallicities against Mg2, a, and environment density 

Next, we investigate whether the mean metallicity of the metal-poor globular clusters correlates 
with other galaxy properties, like metallicity, size or environment. We used the derived mean 
metallicities, as well as the mean (V — I) values (when available), to avoid possible systematic 
errors due to the different transformations from broad-band colors to metallicities. Note that when 
several values for (V — I) were available for a given system, we computed a simple mean, and 
averaged it with any other metallicity determination, if available for that system. 

The correlation between [Fe/H] mp and the Mg2 index of the host galaxy, as well as between 
[Fe/H] mp and the velocity dispersion a of the host galaxy, computed for the full sample, are sig- 
nificant at the >99% confidence level (spearman test returning student distributed t-value= 2.45 
and 2.90 with 35 and 40 degrees of freedom, respectively). Mg2 and a being strongly correlated 
(see Fig. 3 for our sample), the similarity of the relations is not surprising. Neglecting, however, 
the dwarf galaxies (Mg2 < 0.25, a < 150 km-s -1 ) which do not exhibit a separate metal-poor 
component, reduces the significance of the correlations to the <92% confidence level (t = 1.45 and 
1.30 for 33 and 35 degrees of freedom, respectively). Furthermore, the significance is reduced even 
further when selecting the clusters as in Sect. 3.1 (i.e. only galaxies with more than 6 clusters are 
considered). The correlation disappears well below the 90% confidence level (t = 1.07 and 0.75 for 
32 and 34 degrees of freedom, respectively). 

A linear fit returns a relation of the form [Fe/H] mp = 1.07(±0.44) • Mg 2 - 1.68(±0.13) for the 
full sample. The slope of this relation is ~ 10 times shallower than a direct conversion of Mg2 into 
[Fe/H] mp (see Kissler-Patig et al. 1998b) indicating only a very weak dependence of globular cluster 
mean metallicity on galaxy metallicity, if present at all. 

Again, a possible correlation for dwarf galaxies will be discussed in a future paper discussing 
metal-rich globular cluster sub-populations. 

We note that absolutely no correlation is detected with environment density (t = —0.63 for 35 
degrees of freedom). 

3.3. (V — I) against Mg2, er, and environment density 

Similar test as in the above section were preformed for (V — I) instead of metallicity in order 
to avoid any potential systematic effects arising from the conversion of color into metallicity. 
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The most significant correlations of the peak color with a galaxy property are again the cor- 
relation between (V — I) and Mg2, and (V — I) and a for the full sample. The quantities are 
correlated at the ~97% confidence level (t = 1.89 and 1.80 for 27 degrees of freedom respectively). 
These correlations are mostly driven by NGC 4458, the only galaxy with 6 or less clusters for which 
a (V — I) peak color is available. Note, however, that NGC 4458 is a rather uncertain detection: 
Neilsen & Tsvetanov (1999) had 17 clusters to detect a bimodality in the color distribution, of which 
6 are associated to the metal-poor peak. Neither Kundu (1999), nor we could detect a bimodality 
and reproduce this result with 33 clusters detected in the Gebhardt & Kissler-Patig (1999) dataset. 
Removing this galaxy from the sample, the significances of the correlations are reduced to the ~94% 
and ~90% confidence level (t = 1.67 and 1.39 for 26 degrees of freedom, respectively), becoming 
marginal. A linear fit returns a relation of the form (V — I) = 0.24(±0.18) • Mg 2 - 0.89(±0.06) 
for the sample excluding NGC 4458. The slope is compatible with within 1.3 a errors (a similar 
result is obtained for the velocity dispersion). Figure 4 shows the relations. 

A correlation with environment density is completely absent (t = 0.08 for 26 degrees of free- 
dom). 

We conclude that, the mean metallicity of the metal-poor globular clusters is not significantly 
related to the size, metallicity or environment of the host galaxy. A weak correlation might exist 
but remains to be confirmed. 

Interestingly, Larsen et al. (2001 in preparation) find for a very homogeneous sample of 12 
galaxies a similar trend. Tentative implications are discussed below. 

4. Some constraints on the formation of halo globular clusters 

From the results of Sect. 3 we retain two important points : 

1) the mean metallicity of our sample of metal-poor globular cluster systems is only weakly (if 
at all) dependent of the host galaxy properties (My, type, environment, metallicity). This suggests 
that the formation of metal-poor globular clusters was largely uncorrelated with the final host 
galaxy properties. The metal-poor globular clusters (often associated with the halo) could have 
formed either in the proto-galactic phases of the host galaxies or in the earlier phases of the galaxy 
formation. The old age of the globular clusters also supports these latter ideas. 

2) the mean value of our sample of metal-poor globular cluster systems seems almost "uni- 
versal" at [Fe/H] m p ~ — 1.4 with a low dispersion of the metallicity distribution function of 
® [Fe/H] ~ 0.3. This suggests a "universal" mode of formation for the metal-poor globular clusters, 
in the sense that their formation sites had very similar metallicities / properties (?). 
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4.1. Size of the putative fragments 

With this formation hypothesis in mind, we can explore a noteworthy consequence on the 
fragment sizes. Ashman & Bird (1993) already addressed this problem comparing globular clus- 
ter sub-groups within the M31 halo with expectations from cold dark matter models (predicting 
substructures of the order of 10 _3 Mh a i o )- They found the observations and predictions in good 
agreement, with typical mass scales for substructures within M31 of 2 • 10 9 Mq. The sizes can also 
be derived by combining i) the mean metallicity of the metal-poor globular clusters around giant 
galaxies and ii) an assumed relation between the galaxy initial luminosity and the average globular 
cluster metallicity (Cote et al. 1998). Relating the two, most metal-poor globular clusters must 
have formed in objects with luminosities around My ~ —17, i.e. masses of the order of 10 9 to 10 10 
Mq (i.e. in galaxies larger than the remaining dwarfs observed in the Milky Way neighborhood). 
In good agreement with Ashman & Bird (1993) 's results. 

4.2. Was the formation of the metal-poor globular cluster independent of the final 

host galaxy ? 

The current data favor a formation of most metal-poor globular clusters in very similar en- 
vironments/substructures, out of low- metallicity gas (as already speculated by Ashman & Bird 
1993). 

Whether these substructures were fragments a la Searle & Zinn (1978), i.e. entities within 
the dark matter halo of the final host galaxy, or satellites with similar properties but independent 
of the dark matter halo of the final galaxy, is unclear. But it appears secure that they formed 
independently of the metal-rich component (bulge) of the host galaxy. 

Thus, the properties cannot distinguish yet between a scenario in which the metal-poor globular 
clusters formed completely independently of the final galaxy and were later accreted (e.g. Richtler 
1994; Cote et al. 1998; Hilker et al. 1999), and a scenario in which they formed as part of the 
galaxy during the assembly of the halo (Searle & Zinn 1978). The latter would be favored if a 
correlation between the mean metallicity of the metal-poor clusters and the metallicity of the host 
galaxy would exist. The former would be favored if no such correlation would be present. Both 
scenarios would not restrict any formation scenario of the final galaxy (e.g. major collapse, major 
merger, ...), since all galaxy formation scenarios envision similar assemblies of the halos (be it as a 
first stage of further collapse or as the halos of progenitors of a subsequent major merging event). 

If existent, the correlation between the mean metallicity of the metal-poor population and the 
galaxy metallicity is about 10 times shallower than a one to one relation. Such a correlation would 
imply that a fraction of the metal-poor globular clusters formed in satellites (i.e. not related to the 
dark matter halo of the galaxy) but get accreted later while the other part formed in fragments 
(i.e. already within the dark halo of the final galaxy). For our dataset such a correlation is, however, 
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still uncertain, and the fraction of metal-poor clusters from both origins in a system that contributes 
to it is also unknown and probably variable with galaxies. 

Finally, we do know that, in the Milky Way, some halo globular clusters formed around other 
galaxies and were accreted later on. For example, such an accretion process can be witnessed today 
in the form of the Sagittarius dwarf galaxy (Ibata et al. 1994). Although halo globular clusters of 
the Milky Way form a homogeneous population from their metallicity, it has been suggested from 
an analysis of their horizontal branch types (Zinn 1993) that it may contain two sub-systems with 
similar average metallicities. This would complicate the interpretation of their origin. A review of 
pros and cons can be found in Ashman & Zepf (1998), Harris (2000) and Parmentier et al. (2000). 

5. Time and site of formation of the metal-poor globular clusters 

5.1. DLA systems as the progenitors of metal-poor globular clusters 

Pettini et al. (1999) presented a diagram giving the rough location of different components of 
the early universe at z ~ 3 in a N(HI) - metallicity plane. Among the objects whose metallicity can 
be estimated at high redshift are the Damped Lyman a (DLA) systems, the Lyman Break Galaxies 
(LBG) and the Lyman a forest. DLA systems are neutral gas objects observed at all redshifts 
(Pettini et al. 1999; Prochaska & Wolfe 2000). Their dynamics is consistent with protogalactic 
clumps or progenitors of present-day galactic disks (Wolfe et al. 1995; Katz et al. 1996; Haehnelt 
et al. 1996). One interpretation of DLA systems is that they are gas clouds within protogalactic 
halos, that could be associated with Searle-Zinn fragments. LBGs are star-forming objects similar 
to our local starbursting galaxies (Steidel et al. 1996) They appear as objects with compact cores, 
some with multiple components, surrounded by diffuse and asymmetric halos (Steidel et al. 1996; 
Lowenthal et al. 1997). If the metallicity of the Lyman a forest is definitely too low, DLA systems 
and LBGs are found to lie in the same metallicity range as halo globular clusters. However, recent 
data from Kobulnicky & Koo (2000) show that the mean interstellar medium metallicity in their 
LBGs is consistent with the metal-rich globular cluster population in the Milky Way. Consequently, 
since metal-poor globular clusters are formed from metal-poor gas, DLA systems appear as the best 
candidates for their site of formation. 

Furthermore, globular clusters contain stars and the column density of their progenitors must 
be above a threshold of N(HI) ~ 10 20 cm -2 for the star formation to occur (Kennicutt 1989). 
Simulations carried out by Nakasato et al. (2000) also suggest that a self-enrichment in a cloud 
appears to exclude the formation massive star clusters. However, this point is still controversial 
(Parmentier et al. 2000). 

Bringing the above mentioned facts together confirms the allowed location of halo globular 
cluster progenitors in the N(HI) vs. [Fe/H] diagram and suggests that we should concentrate, as a 
working hypothesis, on DLA systems as the potential sites of formation of the (halo) metal-poor 
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globular clusters. 

5.2. Estimating the metallicity evolution of DLA systems 

With the previous hypothesis that DLA systems can be connected to the progenitors of the 
metal-poor globular clusters, we take advantage of the available database of DLA observations, 
normalized to Anders et al. (1989) solar abundances : Log [Fe/H]0 = —4.49 and Log [Zn/H]0 = 
—7.35, to study the chemical evolution of these objects with redshift. However, [Fe/H] may not 
be a reliable estimate of the metallicity of DLA systems, since some Fe may be locked up in dust 
and result in biased [Fe/H] measurements. Pettini et al. (1997) showed that [Zn/H] is a more 
reliable estimator because it essentially measures the metallicity independently of dust depletion. 
The [Zn/H] and [Fe/H] variations as a function of the redshift (Table 3) are plotted in Fig. 5. We 
use the column-density weighted abundances : [< M/Hdla >] = log < (M/H)dla > — log(M/H) Q 
where < (M/H)dla > (M = Fe or M = Zn) and the associated standard deviations as defined in 
Pettini et al. (1997). 

The average values < [Fe/H] DLA >= -1.53 ± 0.40 and < [Zn/H] DLA >= -1.13 ± 0.38 for 
the two complete samples, give < [Zn/ Fe\DLA >= 0.40 ± 0.55 over the whole sample. We do not 
normalize the [Fe/H] to the [Zn/H] values since it remains unclear whether < [Zn/ 'Fe]oLA > is 
constant with redshift. A similar trend of decreasing metallicity with redshift appears both for 
[Fe/Hjoi^ and [Zn/H]^^. The possibility that part of this trend is caused by a fraction of high 
redshifts DLA systems missed because of dust is investigated and ruled out by Pei & Fall (1995) 
(but see also Prochaska & Wolfe 2000). The steeper slope at z > 2.8 has been interpreted as the 
fast formation of metals after the dark age (Pettini et al. 1997; Lu et al. 1996). 

In order to verify whether the above trend of decreasing metallicity with increasing redshift is 
supported by other evidence, we further compute the evolution of the metallicity of the universe 
with the redshift. We use the star formation evolution of (Steidel et al. 1999) and also consider 
an alternative scenario between z = 1 and 3 from sub- mm data (Barger et al. 1999). The second 
scenario implies a higher metal production at z < 3 and a lower metal production at z > 3. The 
expected metallicity is computed following Pettini et al. (1997) except that the values are normalized 
to give a global metallicity in the present-day Universe Z = 1/3 Z (Renzini 1999). As noticed 
by (Pettini et al. 1997), the metal production is deduced from the radiation essentially emitted by 
massive stars. To compare these values with direct DLA metallicity, we need to correct them using 
[a/Fe]= 0.25 (Boesgaard et al. 1999). In our considered metallicity range ([Fe/H] mp < —1.0), this 
value is found approximately constant (Clementini et al. 1999). The resulting curves were added 
in Fig. 5. 

On the one hand, the comparison of DLA and globular cluster system metallicities is direct; 
both are observables. On the other hand, a number of assumptions have been used to estimate the 
metal production in the Universe from the LBGs. Consequently, even if these latter curves may be 
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used as an important check of the DLA metallicity trend, we need to rely on the latter to settle 
our conclusions. 

5.3. A redshift range for the formation of the oldest globular clusters 

The chemical evolution of DLA systems is below the lower limit for our metal-poor globular 
clusters at z ~ 4. The conclusion suggested by these data is that the old globular cluster formation 
occurred at z < 4 for the adopted cosmology. 

Steidel et al. (1999) found that the total integrated UV luminosity at z ~ 3 is of the same 
order as that at z ~ 4 suggesting a similar stellar formation in the two redshift ranges. However, an 
interesting point to stress is that a small number of star-forming galaxies are observed at redshifts z 
> 4 (e.g. Dey et al. 1998; Spinrad et al. 1998; Hu et al. 1999). The observed signal-to-noise ratio of 
these observations is low and we have only a very limited information on the galaxies. However, a 
preliminary conclusion is that those few galaxies may be in the very early stages of their formation. 
Nevertheless, observations of the bulk of high redshift ellipticals are consistent with a formation 
at redshifts of the order of z ~ 3 to 4, while a formation at z < 2 and z > 5 appears to be ruled 
out (Treu & Stiavelli 1999; Menanteau et al. 1999) but see Jimenez et al. (1999) for an alternative. 
To date, the detection of z > 4.5 star-forming galaxies is still anecdotal and most of the stellar 
formation is observed below this redshift. The results presented in this paper brings an additional 
argument to this hypothesis. 

6. Conclusion 

The mean metallicity of metal-poor globular clusters is approximately constant in all galaxies 
and environments, with no significant dependence of galaxy size of metallicity. This argues for 
a formation of all metal-poor globular clusters in very similar gas fragments. Further, it suggest 
either a very homogeneous metallicity of the initial gas out of which old metal-poor globular clusters 
formed, and/or very similar self-enrichment processes of the clouds. Self-enrichment, however, is 
unlikely to play an important role during the formation of clusters, (e.g. Ashman & Zepf 1998; 
Nakasato et al. 2000), so that a very homogeneous metallicity in the initial fragments is favored. 

A weak correlation (to be confirmed) of the mean metallicity of metal-poor globular cluster 
systems with the host galaxy metallicity /size might exist. This would suggest that (at least some 
of) the fragments in which the metal-poor globular clusters formed were already embedded in the 
dark halo of the final galaxy, rather than belonging to independent satellites. 

We found high redshift DLA systems (having high column densities of neutral gas, and similar 
metallicities to the metal-poor clusters and "population II" objects) to be good candidates for the 
formation sites of metal-poor globular clusters. This would support a picture in which, at least 
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some, DLA systems are gas clouds within protogalactic halos. 

We would like to thank Karl Gebhardt for his help in handling the data of the metal poor 
populations and Max Pettini for helpful discussions on DLA systems. Thanks also to Jane Eskdale 
for revising our English in the final manuscript. Last but not least, we thank the referee Keith 
Ashman for very helpful comments that improved the paper. 
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Fig. 1. — Distribution of mean metallicities for the globular cluster systems. The left-hand panel 
presents the distribution of the mean metallicities of metal-poor clusters for the complete sample 
and the right-hand panel only shows globular cluster systems with Nqc > 6 (see text). 
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Fig. 2. — Mean metallicity of the old, metal-poor globular cluster systems plotted against the 
absolute magnitude My of the parent galaxy. Circles represent E/SO galaxies while diamonds 
show spirals. Filled symbols represent globular cluster systems with an observed number of globular 
clusters Nqc > 6, open symbols stand for systems with Nqc < 6. The two lines represent linear 
fits to the complete sample and to the sub-sample with Nqc > 6 respectively. Note that their 
slopes are not significantly different from zero within the errors. 




Fig. 3. — The figure shows the clear correlation between Mg2 and a for the galaxies of our sample, 
as expected Dressier et al. (e.g. 1987). 
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Fig. 4. — The figure shows the peak (V — I) color of the metal-poor globular cluster systems plotted 
against the Mg2 index (left panel) and velocity dispersion (right panel) of their host galaxies. The 
dashed lines show a linear least square fit to the data (excluding the datapoint of NGC 4458). A 
weak (< 2a confidence) correlation might exist between the plotted values. 
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Fig. 5. — Comparative variation of the metallicities with the redshift (Ho = 50km-s _1 -Mpc _1 and 
qo = 0.5). Right hand panel: globular cluster system metallicity distribution. Left hand panel : 
the limits of the globular cluster system mean-metallicity distribution (Nqc > 6) are reported on 
the left hand plot as dashed lines. The uppermost curves have been deduced from Steidel et al. 
(1999, Fig. 9) star formation history (continuous line); the dashed line includes Barger et al. (1999) 
FIR data. [Zn/H] (continuous line crosses) and [Fe/H] (dashed crosses) values of DLA systems 
are taken from Table 3. The age of the oldest Galactic globular clusters is assumed to be > 10 
Gyr, and is reported as an horizontal right-bound arrow. Pettini et al. (1997) noted that assuming 
q = 0.01 would shift the [M/H] by a factor of 2. 
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References. — AB93: Ashman & Bird (1993); BOO: Brodie et al. (2000); B98: Buonanno ct al. (1998); C98: Caldwell et al. (1998); 
D96a: Durcll ct al. (1996a); D96b: Durell et al. (1996b); F96: Forbes et al. (1996); F97a: Forbes et al. (1997a); F97b: Forbes ct al. 
(1997b); F97c: Forbes ct al. (1997c); G96: Geisler et al. (1996); GK99: Gcbhardt & Kissler-Patig (1999); H92: Harris et al. (1992); H00: 
Harris (2000); KP97b: Kissler-Patig ct al. (1997b); K98: Kundu & Whitmorc (1998); K99: Kundu et al. (1999); K99p: Kundu (1999); 
M96: Mighcll ct al. (1996); M98: Montegriffo et al. (1998); N99: Neilsen & Tsvctanov (1999); 098: Olsen et al. (1998); P99: Puzia et al. 
(1999); SGHS98 : Sarajedini et al. (1998); S95: Seeker et al. (1995); S92: Suntzcff et al. (1992); SRN98: Smith ct al. (1998); T95: Testa 
ct al. (1995); W92: Walker (1992); Z95: Zepf et al. (1995) 

Note. — The number of globular clusters used was not recorded by Kundu (1999) but it is typically of the order of 100 or more and 
never below 20. 
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Table 2. Galaxy properties for our sample. Mg2 and a where taken from the HYPERCAT database 
(Prugniel & Simien 1996; Golev & Prugniel 1998), the environment density p from Tully (1988). 
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Table 3. Column Density - Weighted [Zn/H] and [Fe/H]. The [Ln/Yl] DLA are from Pettini et al. 
(1997) while the \Fe/~H]oLA are compiled from Boisse et al. (1998); Lu et al. (1996); Pettini et al. 
(1999); Prochaska & Wolfe (1999). Columns. #5 and #10 give the number of DLA metallicities 
and of corresponding limits for [Fe/H] and [Zn/H] respectively. 



z 


dz 


T 

A mm 


Tmai 


DLA 


[Fc/U] DLA 




[Zn/H] DLA 


V[Zn/H\ 


# DLA 


[Zn/Fe] 


ff [Zn/Fe] 


0.45 


0.1 


5.1 


5.9 


1 (1) 


-1.30 














1.00 


1.0 


5.9 


9.7 


10 (1) 


-1.54 


0.42 


-0.98 


0.33 


4(0) 


0.56 


0.53 


1.75 


0.5 


9.7 


10.5 


5(0) 


-1.32 


0.57 


-0.96 


0.44 


8(2) 


0.36 


0.72 


2.25 


0.5 


10.5 


11.0 


13 (0) 


-1.61 


0.50 


-1.23 


0.38 


12(6) 


0.38 


0.63 


2.75 


0.5 


11.0 


11.4 


4(0) 


-1.29 


0.64 


-1.11 


0.27 


7(4) 


0.18 


0.69 


3.25 


0.5 


11.4 


11.6 


3 (0) 


-1.75 


0.10 


-1.39 




3(3) 






3.75 


0.5 


11.6 


11.8 


2 (1) 


-1.95 


0.20 












4.25 


0.5 


11.8 


12.0 


3(1) 


-2.28 


0.25 












2.45 


4.1 


5.1 


12.0 


47 (4) 


-1.53 


0.40 


-1.13 


0.38 


34(15) 


0.40 


0.55 



